Gravitational collapse of dustlike matter with heat flux 
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We present a new class of solutions to Einstein equations for the spherical collapse of dustlike 
matter coupled with heat flux. In this family of solutions spacetime shear is necessarily non-zero. 
Also these solutions have an interesting property that there is always a bounce before the singularity, 
which is caused entirely due to the dissipative processes. We show there exist open sets of initial 
data for which the bounce occurs before any trapped surface formation, making the star explode 
away to infinity. We also discuss the role of heat flow in generating spacetime shear and in modifying 
the effective inertial mass of the matter cloud. 
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I. INTRODUCTION 

A considerable amount of work has continued in re- 
cent years to model a realistic gravitational collapse sce- 
nario within the framework of Einstein's gravity (see e.g. 
[H, and the references therein for some recent re- 
views). The reason for this interest is that general rela- 
tivity generically admits the existence of spacetime sin- 
gularities. These are extreme regions in the spacetime 
where densities and spacetime curvatures typically blow 
up and the theory must breakdown. 

When a sufficiently massive star starts collapsing grav- 
itationally on exhausting its nuclear fuel, it would not 
settle to a stable configuration such as a neutron star. 
What happens in such a case is an continual gravita- 
tional collapse ensues, where the sole governing force is 
gravity. One of the interesting problems in this pursuit of 
study of gravitational collapse of such a massive star, is 
studying various dynamical solutions of Einstein's equa- 
tions with radial heat flow in the interior, satisfying all 
the energy conditions. This radial heat flux is generated 
by various dissipative processes in the collapsing interior 
and indeed the relevance of these processes are central to 
a realistic astrophysical collapse. 

The earlier investigations done in this context, sug- 
gest that the final outcome of continual collapse in the 
presence of a radially outgoing heat flux is dramatically 
different from the usual Oppcnheimcr-Snyder picture in 
which a trapped surface develops before the singularity 
and ultimately the collapsing matter settle down to a 
black hole. For example, the solution presented in Q 
describes the collapsing matter which radiates out mass 
due to the outward heat flow to an exterior Vaidya re- 
ion (for details of the matching conditions please see 
) and through this process the complete star evapo- 
rates away leaving behind a flat Minkowski spacetime. A 
different solution was found in 1 and @, where the loss 
of matter due to heat flux stops the trapped surface for- 
mation and a naked singularity is formed as the end state 



of the gravitational collapse. Detailed investigations on 
the dissipative processes in the interior spacetime, using 
diffusion approximation was carried out in Q- [ill ], and 
it was shown that due to the inertia of heat (as described 
by Tolman [12]); the collapse may slow down and stop 
before the singularity formation. 

However, in all these earlier works we referred to above, 
the interior spacetime was considered to be shearfree, 
which is rather a strong supposition. The basic aim of 
this paper is to investigate the final outcome for a more 
general class of collapsing solutions of Einstein equations 
for a spherically symmetric spacetime, without the re- 
striction of being shearfree. For better clarity we take 
the matter to be dustlike (i.e the pressures are small 
compared to the density at all epochs) coupled to a ra- 
dially outgoing heat flux and satisfying all the energy 
conditions. We discuss various regularity issues and do a 
near the center analysis of the class of solutions presented 
here, to show that in presence of heat flow the collapsing 
configuration always undergo a 'bounce' before the sin- 
gularity and hence the singularity formation is avoided. 
We also discuss the role of heat flow in generating space- 
time shear and in modifying the effective inertial mass of 
the matter cloud. 



II. THE COLLAPSING SYSTEM 

A general spherically symmetric metric in the coordi- 
nates (t, r, 9, <fi) can be written as, 

ds 2 = g ab dx a dx b + R 2 (t, r)dfl 2 (1) 

where where a,b run from to 1, x° — t, x 1 = r and 
dfl 2 is the line element on a two-sphere. In general, we 
have freedom of coordinate transformations of the form 
t' = f(t,r) and r' = g(t,r). We use one transformation 
to make the gt r term vanish. So the metric now take the 
form, 



ds 2 



(2) 



'Electronic address: goswami@phys.ualberta.ca 



The energy-momentum tensor for the dustlike matter 
coupled with heat flow is given as 

rpik 



pu i u k 



q l u k 



u i q k 



(3) 



2 



where p is the energy density, u 1 is the timelike unit vec- 
tor (or the 4- velocity of the collapsing fluid) and q l is the 
heat flow vector. We now use the other transformation 
to make the 3-velocity of the collapsing dust to vanish 
{i.e. the 4- velocity of the collapsing dust is of the form 
u l = e~"<5|), so that the frame is comoving Now we 
do not have any gauge freedom of two variables left. We 
still have the two scaling freedoms of one variable, which 
we would use later. It is easily seen from (|3j) that by re- 
defining the density we can always make the heat flow to 
be orthogonal to the velocity vector such that q l Ui = 0. 
In that case the vector q 1 would have the form, 



q l = Q(t,r)5l 



(4) 



The matter cloud is assumed to have a compact support 
with < r < rf,, where rj denotes the boundary of the 
cloud, outside which it is to be matched through suitable 
junction conditions with an exterior spacetime geometry. 
For example, as we would like to model an astrophysical 
collapse scenario here, we would match the interior with 
an outgoing (exploding) generalized Vaidya exterior. The 
function R(t, r) is the area radius of a shell labeled V at 
an epoch 't'. And as we consider the situation where 
the cloud is collapsing initially we must have R < at 
the initial epoch. In a continual collapse scenario this 
ultimately reaches R = 0, which is the singularity where 
all matter shells collapse to a zero physical radius. We use 
the scaling freedom for the radial coordinate r to write 
R = r at the initial epoch t — U, and to differentiate 
the center of the cloud from the singularity (as in both 
these cases the area radius vanishes) introduce a function 
v(t,r) as defined by 



v(t, r) = R/r . 



(5) 



We then have 



R(t, r) = rv(t, r), v(U,r) = 1, v(t s (r), r) = (6) 

with v < 0. The time t = t s (r) corresponds to the shell- 
focusing singularity at R = 0, where the matter shell 
labeled V collapse to a vanishing physical radius. We 
assume that the matter field satisfies all the energy con- 
ditions, in terms of the density and heat flow vector this 
corresponds to [4| 

P>2\q\ (7) 

where \q\ — gikq l q k ■ The dynamic evolution of the initial 
data is determined by the Einstein equations. For the 
metric ([2]), and with the following definitions (to make 
the second order system of equations to a coupled first 
order system), 

G(t, r) = e- 2 ^(R'f, H(t, r) = e- 2v (R) 2 , (8) 



F = R(l-G + H), (9) 



we can write the independent Einstein's equations for the 
given matter field ^ as, 



F' = pR 2 R', F = , 



v R p = — e 



-^ +u) d t [e^QR 2 ] , 



R'% - 2Rv' = Re^+^Q . 



(10) 



(11) 



(12) 



Here (') denotes the partial derivative with respect to the 
co-ordinate V while (') denotes partial derivative with 
respect to the co-ordinate 't\ The function F here has 
an interpretation of the mass function for the cloud, and 
it gives the total mass in a shell of comoving radius V 
on any spacelike slice H = const' '. The energy conditions 
imply F > 0. To preserve regularity at any non-singular 
epoch 'f, we must have the mass function should vanish 
at the center of the cloud. 

As in Lemaitre Tolman Bondi (LTB) models, we now 
try to find a new family of solutions with v = v{t). This 
is an ansatz that we use here and this would only be 
consistent with some classes of heat flow and not the 
most general class. Then we can always rescale the time 
co-ordinate to make the metric component gu = — 1. We 
note that then we have no more gauge or scaling freedom 
left. Now from ITU1) we have, 



(13) 



F = F(r) = r 3 M(r) 



where M. (r) is at least a C 2 function of the co-ordinate 
V. It is to be noted that F must have this form, which 
follows from the regularity and finiteness of the density 
profile at the initial epoch t = ti, and at all other later 
regular epochs before the cloud collapses to the final sin- 
gularity at R = 0. Also integrating (JTTJ) we have, 



Q = 



5 /(r)e 



-3V 



R 2 



(14) 



where f(r) is another C 2 function, and the form of the 
function of integration is chosen in a way such that the 
heat flow vanishes at the center of the cloud. We again 
would like to emphasize that this is not the most general 
class of the heat flux function. However as we would see 
later, near the center of the cloud and at any given epoch 
v, Q « rp. This definitely makes sense physically, as the 
heat flux must vanish at the center and also it should 
be increasing with increase in density. Substituting the 
form of v and (Tj"4"| in (TT2"]) , we have 



R G- 



R 



(15) 



Multiplying both sides of the above equation by 
— (l/2)(l/i?' 2 ) and rearranging we get 



9t[G ( 



-1/2)1 



1 r 3 /(r) 
2R{R') 2 



(16) 
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Now to construct classes of solutions to Einstein equa- 
tions, let us define a suitably diffcrcntiable function 
A(r,v) as follows, 



1 r 3 f(r) 
2R{R') 2 



A(r, v) 



(17) 



We would like to emphasize here that this is no special 
choice but merely a construction which would help us 
decouple the Einstein equations and find classes of solu- 
tions. Using this definition we can now integrate (flu]) to 
get 



G 



1 



[A(r, v) + r 2 b(r)} 2 



(18) 



Here again, the function b(r) is a suitably differentiable 
function of integration. We also would like to point out 
here that the function A(r, v) = A(r) corresponds to 
f(r) = (that is the situation where Q = 0) and we 
get back the well known LTB solution. Now from © we 
see that in order to make v well defined at the center of 
the cloud (r = 0), we should have [A(r, v) + r 2 b(r)](~ 2 ' 
1 + r ' B(r,v), where B(r,v) is another suitably differen- 
tiable function. Or in other words the function A(r, v) 
must have the form, 



A(r,v) 



1 



y/1 + r 2 B(r,v) 



2 b(r) 



(19) 



We note that the above equation is not really any ansatz 
or a special choice, but quite a generic class of function 
consistent with and allowed by the required regularity 
conditions of metric functions at the center of the cloud. 
Comparing (fl8)) with flSJ and using (fl~9| we get 



_ „2ifi _ 

Qrr 6 



(R)' 2 



1 + r 2 B(r,v) 
Using the above equation and © we have, 



B(r,v) 



M{r) 



(20) 



(21) 



The negative sign in the above equations corresponds to 
a initial collapse scenario, where we have R < initially. 
Also using equation (fT!)]) in (TTT)) we get, 



v + rv 



f(r)(l + r 2 B(r, v))( 3 / 2 ) 
vB(r,v), v ^B(r,v) + ^ 



(1/2) 



(22) 



Now from equations (|21| and (|22|) we get the expressions 
of v' and v in terms of r, v and the function J3(r, v) and 
it's partial derivative with respect to v. Thus, to get a so- 
lution of v(t, r) we need the Pfaffian differential equation, 
v'dr + vdt — dv — , to be integrable. This integrability 
conditions gives us the required second order differential 
equation for the function B{r, v). To show that solutions 



of this equation exist, it is sufficient to find one solu- 
tion. We have already seen that B = B(r) is a solution 
which corresponds to the LTB solution with f(r) = 0. In 
that case, providing the value B(r,l), completely deter- 
mine the allowed function B(r,v). From ([22]) it is clear 
that we get B(r, v), v \ v= \ (which corresponds to providing 
Q.t(U, r)) from the other three free functions. 

Thus, we see that in this class of models there are 
three free functions of the co-ordinate V, namely M.(r), 
B(r,l) and f(r). The first function gives the density 
profile, the second one describes the velocity profile of the 
infalling shells while the third describes the flux profile at 
the initial epoch. Given these initial profiles we can solve 
the system of Einstein equations completely, and from 
pTjl we can write the solution of v(t,r) in the integral 
form as, 



t(v,r) = 



y/vdv 



v y/vB(r,v) +M(r) 



(23) 



We note that the variable V is treated as a constant in 
the above equation. The above equation gives the time 
taken for a shell labeled V to reach a particular epoch 
from the initial epoch v = 1. 

To complete a collapse model, we need to match the 
interior spacetime to a suitable exterior spacetime. As 
we are interested here in modeling collapse of astrophys- 
ical objects (such as massive stars), we have assumed the 
matter to have compact support at the initial surface, 
with the boundary of the cloud being at some r = r^. 
We match the interior with a general class of exterior 
metrics, which are the generalized Vaidya spacetimes (it 
describes the matter which is a combination of matter 
fields of Type / and Type II [3])- at the boundary hyper- 
surface S given by r = rb. For the required matching we 
use the Israel-Darmois conditions ([H[- p7|). where we 
match the first and second fundamental forms (the met- 
ric coefficients and the extrinsic curvature respectively), 
at the boundary of the cloud. The metric in the exterior 
of T, is given by, 

ds 2 + = -(l- 2M ^ V A dV 2 -2dVdr v +r 2 v dQ 2 (24) 



where V is the retarded (exploding) null co-ordinate and 
r v is the Vaidya radius. For details of matching condi- 
tions we refer to [21 . 



III. ANALYSIS OF THE SOLUTION NEAR THE 
CENTER 

To understand the dynamics of the collapsing system 
in more transparent way, let us consider a small neighbor- 
hood of the center of the cloud. That is we consider the 
value of the coordinate V to be small enough so that we 
can neglect the expressions containing the higher pow- 
ers of r. Let us consider the free functions having the 
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FIG. 1: The function Bq(v) as found by numerical integration, 
for different initial values. The constants were taken as |/o| = 
0.5, Mo = 1 We can see that even if we start from a positive 
value at the initial epoch, the function becomes negative at 
the later epoch causing a bounce. 



FIG. 2: The plot of the vv 2 as a function of v for different 
initial values. The constants were taken as |/o| = 0.5, Mo = 1. 
It is seen that before the spacetime singularity v = 0, the 
collapsing velocity goes to zero causing a bounce. 



following Taylor expanded form near the center, 

f(r) = fo + / 2 r 2 + • ■ • 
M(r) = M + M 2 r 2 + ••• (25) 

We also write the function B(r, v) in the form, 

B(r,v) = B Q {v) + B 2 (v)r 2 + ••• (26) 

We know that near the center R' w v. Putting the above 
equations in ([22)) and neglecting higher powers of r, we 
get, 



(27) 



Now since we consider the heat flux is outwardly directed 
we must have /(r) < for all r and hence fo < 
Hence we get the following ordinary differential equation 
for the evolution function Bq(v), 



B (v) jV 



l/o 



vW>y/B (v)v + M 



(28) 



It is not possible to integrate the above equation an- 
alytically in a closed form. However numerical integra- 
tion shows that even if we choose Bq(1) > at the ini- 
tial epoch, the function Bq(v) changes sign at some later 
epoch and there always exist an epoch Vb € (0, 1] where 
V = 0. In other words the collapse of the small neighbor- 
hood of the central shell stops and shells bounces back 
before the spacetime singularity. Figure 1 shows the func- 
tion Bq(v) as a function of v and we see that there is 
an epoch at which Bq(v), v blows up. From the above 
equation it is clear that it is the same epoch in which 
\J 'B {v)v + M — ► or the quantity v — (as shown in 
Figure 2) and the shells rebounce. Though this analysis 
is for a small neighborhood around the center, we note 
from equation (|22p that we are only considering those so- 
lutions where R' > throughout the spacetime. In other 



words we are considering solutions free of any shell cross- 
ing singularities. In this scenario, if there is a bounce 
around the center, we can conclude the the whole col- 
lapsing matter bounces back and there are no spacetime 
singularities. 

From equation 1(7} we see that, as \q\ « r 2 near the 
center, all energy conditions are satisfied in the neigh- 
borhood of the center. Further to this, we can always 
choose the free functions F(r), f(r) and B(r, 1) in such 
a way that energy conditions are obeyed throughout the 
spacetime. To inspect more carefully the reason for this 
bounce, in spite of satisfying all energy conditions, let us 
consider the 'acceleration 1 term, 



l/o 



Mo 

2 



u( 5 / 2 ) y/B (v)v + M v 



(29) 



From the above equation it is clear that even if we 
choose the initial data in such a way that at the initial 
epoch v < (i.e. the collapse is accelerating), as the 
collapse proceeds the first term in the equation starts 
dominating and ultimately there is a decelerating effect 
as v > 0. At the epoch of the bounce we have v — > oo, 
that is, there is an infinite deceleration. 

There is also an alternative thermodynamical explana- 
tion to this phenomenon. As shown in [TlT |. if we cou- 
ple the Einstein equations with Israel Stewart transport 
equation (l8j we can easily see that the Effective Inertial 
Mass Density of the infalling matter (as defined by the 
'Newtonian' form of the equation of motion) reduces due 
to the presence of heat flow and becomes negative as the 
collapse proceeds. For the dustlike matter field we con- 
sider here, the effective inertial mass density behaves like 
p(l — a), where the function a depends on the tempera- 
ture, density, thermal conductivity and relaxation time. 
During the collapse the temperature of the cloud may in- 
crease enormously which in turn can make a > 1. Thus 
in spite of the fact that infalling matter should gain mo- 
mentum due to the outgoing flux, the negative effective 
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inertial mass density actually creates a 'force' against the 
motion and the collapse stops before the singularity. 

Coupling of heat flux with collapsing matter also has 
another important physical effect. We know, for the 
given metric ([2]) , the shear tensor for the collapsing mat- 
ter is given by, 

^ = -r-k"(i-*) < 30 » 

Now, in the absence of any heat flow, if we consider 
B(r,v) — 1 and M(r) = Mo, then the spacetime is 
shearfree and the collapse becomes Oppenheimer Snyder 
homogeneous dust collapse. However if we couple the 
matter with heat flux then we see that the shear tensor 
is necessarily non-zero. In other words even if we start 
with homogeneous matter and couple it with heat flux 
then the matter becomes inhomogeneous at later epochs 
and spacetime shear develops. We know that presence 
of shear delays the formation of trapped surface 
Also the loss of matter due to continuous radiation in 
the generalized Vaidya exterior further helps in delaying 
the trapped surface formation. In that case we can easily 
see that there exist a large class of collapsing solutions 
(depending on choice of initial data and) for which the 
collapsing star bounces back before any trapped region 
forms, that is, the star explodes away to infinity. 

IV. CONCLUDING REMARKS 

In this paper we presented a new family of solu- 
tions which generalizes the well known Lemaitre-Tolman- 
Bondi solutions, for the collapse of dustlike matter with 
heatfiux and the spacetime shear is necessarily non-zero. 
The most important feature which emerges from these 



class of solutions is that, there is always a bounce be- 
fore the singularity, which is caused entirely due to the 
dissipative processes. Also we showed, there exist open 
sets of initial data for which the bounce occurs before 
any trapped surface formation, making the star explode 
away to infinity. 

Such a study of gravitational collapse is warranted for 
several reasons. First of all, in the later stages of a mas- 
sive star evolution, due to considerable decrease in opac- 
ity of fluids (for photons and neutrinos), heat conduction 
does play a very important role and that cannot pos- 
sibly be neglected. It is shown, for example, that this 
phenomenon of decreasing opacity give rise to neutrino 
radiative heat conduction [201 ] in the stellar fluid which 
may increase the temperature of the collapsing star to 
a large extent. Secondly, as we have already seen, the 
general properties of the collapse in presence of heat flux 
is dramatically different from when it is absent and we 
should give it a careful consideration, because there may 
exist several alternate pathways through which a mas- 
sive star may actually die. And finally as we know there 
are deep paradoxes associated with spacetime singulari- 
ties where the classical theory breaks down. While one 
expects a successful quantum picture of spacetime would 
actually be able to resolve these singularities, it is always 
welcome if they are removed at the classical level. 
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